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ABSTRACT 

We present a new statistical method to determine the relationship between the 
stellar masses of galaxies and the masses of their host dark matter haloes over the en- 
tire cosmic history from z ~ 4 to the present. This multi-epoch abundance matching 
(MEAM) model self-consistently takes into account that satellite galaxies first become 
satellites at times earlier than they are observed. We employ a redshift-dependent pa- 
rameterization of the stellar-to-halo mass relation to populate haloes and subhaloes 
in the Millennium simulations with galaxies, requiring that the observed stellar mass 
functions at several redshifts be reproduced simultaneously. We show that physically 
meaningful growth of massive galaxies is consistent with these data only if obser- 
vational mass errors are taken into account. Using merger trees extracted from the 
dark matter simulations in combination with MEAM, we predict the average assembly 
histories of galaxies, separating into star formation within the galaxies (in-situ) and 
accretion of stars (ex-situ). Our main results are: The peak star formation efficiency 
decreases with redshift from 23 per cent at z = to 9 per cent at z = 4 while the 
corresponding halo mass increases from lO n - 8 M to 1O 12 5 M . The star formation 
rate of central galaxies peaks at a redshift which depends on halo mass; for massive 
haloes this peak is at early cosmic times while for low-mass galaxies the peak has 
not been reached yet. In haloes similar to that of the Milky- Way about half of the 
central stellar mass is assembled after z = 0.7. In low-mass haloes, the accretion of 
satellites contributes little to the assembly of their central galaxies, while in massive 
haloes more than half of the central stellar mass is formed ex-situ with significant 
accretion of satellites at z < 2. We find that our method implies a cosmic star forma- 
tion history and an evolution of specific star formation rates which are consistent with 
those inferred directly. We present convenient fitting functions for stellar masses, star 
formation rates, and accretion rates as functions of halo mass and redshift. 

Key words: cosmology: dark matter, theory - galaxies: evolution, high- redshift, 
statistics, stellar content 



1 INTRODUCTION 

Modern large-scale dark matter simulations are able to 
make accurate and definite predictions for the prop- 
erties of dark matter haloes at all cosmic epochs in 



cold dark matter (CDM) cosmo lo gies (ISpringel et al] 
20051: iBovlan-Kolchin et alj 120081 : iKlypin et alj l201ll : 



Angulo et alj|2012h . In this framework, galaxy formation is 
driven by the growth of the underlying large-scale structure 
and the formation of dark matter haloes: galaxies form by 
the cooling and condensation of gas at the centres of the 
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poten t ial wells of such dark m atter haloes (|White fc Reesl 
Il97g| : iFall fc Efstathioul |l980T ). This results in a tight 
correlation between the physical properties of galaxies, such 
as their stellar mass and star formation rate (SFR), and 
the mass of the haloes by which they are hosted. 

Galaxy formation studies have adopted several ap- 
proaches to link galaxy properties to those of the dark 
matter haloes. The most direct is to model the phys- 
ical mechanisms that drive the formation of galaxies 
from first principles. Here, an JV-body treatment of the 
dark matter evolution is combined with a treatment 
of the bar yonic component based either on a hydro- 
dynamical (|Katz et alj 1 19961 : ISpringel fc HernquistJ 120031 : 
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Keres et al l 120051; ICrain et all 120091 ; ISchave et all |2010t 

Qppenheimer et al] 2010f) oron a sem i-ana lytic method 



llKauffmann et al.l 119991: ISpringel et al.) 12001 1: iHatton et al 
20031 : ISpringel et al.l 1120051: iKang et all bood; broton et a! 



20061 : iBower et al.ll2006l : iGuo et al.ll201ll )~espite substan- 



tial progress, simplified and uncertain recipes must still be 
employed to model the formation of stars and black holes 
and the associated feedback processes, and the detailed 
properties and evolution of observed galaxy populations are 
still only partially reproduced. 

An alternative approach to link galaxies to haloes 
has emerged with t he accumulation of data from 
large galaxy sur v eys (|Cole et al. _ 200 ll; Bell et ahl 



Drory ct al. 



c y sur v eys (lUole et alj I^UUil: men et alj 
Jl2005l: iFbiitana et alj 120061: iPanter et all 



Perez-G onzalez ct al. 2008; Li & White 120091 : Santini et al.l 



200 



200 



20121 ) . Instead of attempting to model the baryonic physics 



that drive galaxy formation, dark matter haloes are popu- 
lated with galaxies using a simple empirical method with ad- 
justable parameters. In the popular halo occupation distri- 
bution (HOD) formalism, the distribution of galaxies (hav- 
ing certain intrinsic properties such as luminosity, color, or 
type) within main haloes of a given mass is constrained 
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). A variant of the classical HOD approach is the con- 



120031 : Ivan den Bosch et all 120031 : lYang et al] l201ll ). Since 
reliable galaxy clustering measurements are unavailable at 
high redshift, the HOD and CLF approaches have typically 
been used at low redshift. 

In order to circumvent this problem, galaxies can be 
linked to the underlying substructure directly by assum- 
ing a monotonic relation between stellar mass and halo 
mass, and matching the cumulative abundance of galax- 
ies to those o f haloe s and subhaloes Jyalc & Ostrike r]l2004l: 
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The only observational input to this 'abundance 
matching' method is thus the observed stellar mass func- 
tion (SMF) or luminosity function. This approach then pre- 
dicts clustering statistics remarkably well. In this frame- 
work, the stellar-to-halo-mass (S HM) relation can b e ob- 
tained in two ways: The first (e.g. IConrov et al.|[2006T ) is to 
create a mock galaxy catalogue from the observed SMF for 
the simulated volume, to order both galaxies and haloes by 
mass, and then to match them one by one. A fitting function 
can then be used to get the SHM relation. The second way 
(e.g. iMoster et aLll2010h is first to assume a functional form 
for the SHM relation, and then to populate the haloes with 
galaxies and compute the model SMF. The parameters of 
the SHM function are then adjusted so that the observed 
SMF is reproduced. The latter approach has the advantage 
that scatter in stellar mass at a given halo mass can easily 
be added. 

The underlying assumptions in this model are that 
satellite galaxies are living in subhaloes and were central 
galaxies before their own haloes were accreted and became 
subhaloes. The properties of the galaxy population today 
are thus a consequence of the formation of central galaxies in 
main haloes at different redshifts and the subsequent accre- 



tion and evolution of satellites. Since the subhalo is affected 
much more by tidal stripping than the galaxy at its centre, 
the stellar mass of satellites is more closely related to subhalo 
mass at the time of accretion rather than to current sub- 
halo mass. Most abundance matching studies si mply assign 
the in fall mass to subhaloes. As pointed out bv lYang et al.l 
l|201ll ), this is equivalent to assuming that the SHM relation 
does not evolve with redshift, i.e. that the relation between 
the stellar mass of central galaxies and their halo mass to- 
day is the same as the relation between the stellar mass of 
satellite galaxies and their halo mass at the time of accre- 
tion. However, applying the technique a t different redshifts 
results in an evolvi ng SHM relation (e.g. IConrov et al]|2006l ; 
IMoster et al]|201ol ). This indicates that the classical abun- 
dance matching approach is not self-consistent. 

This problem can be solved if we assign the stellar mass 
of the satellites depending on the infall mass of their sub- 
halo and the infall redshift. The stellar mass of the central 
galaxies still depends on the main halo mass at the redshift 
of observation in the same way. This is equivalent to as- 
suming that the stellar mass of a satellite does not change, 
once it has entered the main halo. In order to implement 
this scheme, we first use the classical abundance matching 
approach to see how the parameters of the SHM relation 
depend on redshift and find a redshift-dependent parame- 
terization of the SHM relation m*(Afh,z). We then employ 
this parameterization to put a galaxy at the centre of every 
halo and subhalo in the simulation at all redshifts between 
z = and z = 4, adjusting the parameters of the relation 
to reproduce simultaneously all the observed SMFs in this 
redshift range. In this way we use a self-consistent relation 
between stellar mass, halo mass and redshift for both cen- 
trals and satellites. 

Statistical models, such as the abundance matching ap- 
proach, are typically employed to derive constraints on the 
relation between galaxy properties and halo mass at a given 
epoch. With the advent of wide high-redshift surveys, it has 
become possible to use this connection, together with in- 
formation on the evolution of dark matter haloes, to study 
the e volution of galaxy properties througho ut their history 
(e.g. IWhite et alj|2007l; IConrov et al] 120071 1. In particular, 
IConrov fc Wechslerl (|200gF show how this method can be 
extended to put constraints on the SFR of galaxies. For this 
they use an analytic prescription for galaxy and halo abun- 
dances and the evolution of haloes, focusing on z < 1. We 
follow a similar approach, but use iV-body simulations to get 
halo abundances and assembly histories. Then we identify 
the progenitors of a given halo at an earlier epoch and derive 
the amount of stellar mass that has been accreted onto the 
central galaxy and the total stellar mass by which it grew. 
Finally, employing a recipe to account for stellar mass loss, 
we compute the stellar mass that has been formed through 
star formation and convert this into a SFR. We are thus able 
to derive the star formation history (SFH) of the galaxy at 
the centre of the halo. 

This paper is organized as follows. In Section 2, we de- 
scribe the iV-body simulations including the halo finding 
algorithm and the treatment of 'orphan' galaxies. We also 
summarize the observational SMFs that have been used. In 
Section 3 we describe the details of the model that matches 
galaxies to haloes. In Section 4 we use this relation and the 
halo accretion histories to derive SFR histories as a function 
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of halo mass and we show that our model can reproduce the 
observed evolution of the cosmic SFR density and of the spe- 
cific SFRs of galaxies. We discuss some of the implications 
of our model in Section 5 and summarize in Section 6. 

Throughout this paper, we assume a 
' WMAPT ACDM cosmology with (fi m , fi A , fi b , h , n > °"s) = 
(0.272, 7 28, 0.0 46, 0.704, 0.967, 0.810). We employ a 
IChabrierl (|2003l ) initial mass function (IMF) and we con- 
vert all stellar masses to this IMF. All virial masses are 
computed with respect to 200 times the critical density. In 
order to simplify the notation, we will use the capital M to 
denote dark matter halo masses and the lower case m to 
denote galaxy stellar masses. 



2 SIMULATIONS AND OBSERVATIONS 

In this paper we aim to link observable properties of galax- 
ies to the underlying dark matter distribution using a sta- 
tistical sample of both galaxies and dark matter haloes. The 
properties of the halo population originate fro m two TV-body 
simu lations: the Millennium Simulation (MS ISpringel et all 
120051 ) for massive haloes and the Millennium-II Simulation 
(MS-II iBovlan-Kolchin et al.ll2009l ) for lower mass haloes. 

For the galaxy population a t low redshi f t we use the 
precise stellar mass functions of iLi fc White! (2009) based 
on the spec troscopic data se t of the Sloan Digital Sky Sur- 
vey (SPSS lYork et al.ll2000l ). At high redshift we use two 
sets of o bserved SMFs: for massive galaxies we use the 
SMFs by iPerez-Gonzalez et alj (|2008l . PG08) taken from a 
wide survey, a nd fo r low-mass galaxies we use the SMFs by 
ISantini et alJ ()2012l . S12) taken from a deep survey. 



2.1 Numerical Simulations 

Both Millennium simulations adopt the concordance ACDM 
cosmology with parameters chosen to agree with the first- 
year Wilkinson Mic rowave Anisotropy Probe ( WMAP) data 
ijS pergel et aUl20 03) . As these parameters are not consistent 
with the latest WMAP data an alysis, we use the scaling 
recipe bv lAngulo fc White! (|2010h to convert the properties 
of the simulations to represent the growth of structure in 
a WMAP1 cosmology. Specifically, we chose fl m = 0.272, 
Q.K = 0.728, fi b = .0455, Hp = 70.4 km s" 1 Mpc" 1 , n = 
0.967 and cr 8 = 0.81 (jKomatsu et al.ll201ll ). 

Both simulations were carried out in a periodic box ad- 
vancing 2160 3 dark matter particles from redshift z ~ 127 
to 0. The side length of the box is 761 Mpc for the MS, and 
152 Mpc for the MS-II. The corresponding particle masses 
are 1.58 x 10 9 M Q and 1.27 x 10 7 M Q , respectively. The large 
volume of the MS enables us to sample reliably even the cen- 
tral galaxies of rare and massive clusters, while the MS-II 
has sufficient resolution to sample the faintest dwarf galax- 
ies. 

Dark matter haloes are identified at each output using 
a friends-of-friends (FoF) halo finder that links together par- 
ticles with separation less than 0.2 of the mean interparticle 
distance. Substructur es inside the FoF grou ps are then iden- 
tified using SUBFIND (jSpringel et alJlioOll ). The most mas- 
sive self-bound subgroup is defined as the 'main halo', and 
all other subgroups within the FoF groups are referred to as 
'subhaloes'. The virial radius and mass of the main haloes 



are determined with a spherical overdensity criterion: the 
density inside a sphere centered on the most bound particle 
is required to be equal to 200 times the critical value. Merger 
trees were then built to link each subhalo found at a given 
snapshot (including the main haloes) to a unique descendant 
in the following snapshot. With these merger trees we can 
track the evolution of each halo throughout the simulation. 

In this work, we assume that both main haloes and 
subhaloes host galaxies at their centres. We will refer to the 
galaxy at the center of a distinct halo as a 'central galaxy', 
and the galaxies within subhaloes as 'satellites'. There is 
thus one central galaxy within each FoF group but there 
may be any number of sattelites (including zero). Ab initio 
models of galaxy formation predict that the stellar mass of 
a galaxy is tightly correlated with the depth of the potential 
well of the halo in which it formed. For central galaxies, 
the relevant mass is the virial mass of the main halo at 
the time of observation, but for satellites the situation is 
more complex: Subhaloes lose mass while orbiting in a larger 
system as their outer regions are tidally stripped. Stars are 
centrally concentrated and more tightly bound than the dark 
matter, however, so the stellar mass of a satellite galaxy 
changes only slightly after it enters a larger halo. Therefore 
the subhalo mass at the time of observation is not a good 
tracer for the stellar mass of a satellite. A better tracer is the 
subhalo mass at the time that it was accreted by the main 
halo, which typically corresponds to the maximum mass over 
its history. 

Due to the finite mass resolution of the simulations, sub- 
haloes can no longer be identified once their mass has been 
tidally stripped below the resolution limit. Since mass loss 
can be substantial this is important even for fairly massive 
subhaloes. A special treatment of these so-called 'orphans' is 
necessary. We determine the orbital parameters the last time 
a subhalo is identified in the s imulation and use them in the 
dynamical friction formula of iBinnev fc Tremaine! |l987): 



GA'/ sa t In A 



(1) 



where r sat is the distance between the centers of the main 
halo and of the subhalo at the time the merger clock was 
set, M sat is the mass of the subhalo at this time, and 
In A = ln(l + Myjr/Msat) i s the Coulomb logarithm. Fol- 
lowing |Pe Lucia fc Blaizotl (|2007T ) we set ayf = 2.34, which 
was also found to be a pprop riate in the A^-body studies by 
iBovl an-Kolchin et alJ (|2008l ). We retain the disrupted sub- 
halo for the purposes of abundance matching until a time t<n 
has elapsed, and assume that only then does it truly merge 
with the main halo. 



2.2 Observed stellar mass functions 

Having derived the abundance of the halo/subhalo popu- 
lation as a function of mass and redshift, we now have to 
select samples of observed galaxy abundances, to link the 
two populations. At low redshift (z = 0.1), the SDSS galaxy 
catalogues allow a determination of the local SMF down 
to a stellar mass of 1O 7 M0. Above about 1O 8 M0 it is ro- 
bust against incompleteness. The large volume of the sam- 
ple leads to very low uncertainties due to cosmic variance 
jTrenti fc Stiavellill2008l ; iMoster et al.ll201ll ). 

We employ the recent measurements by ILi fc White! 
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(2009), base d on a complete and u niform sample from the 
SDSS/DR7 jAbazaiian et al.| [2009V and we use the conver- 
sion to stellar mass ba sed on the S PSS r-band model lumi- 
nosities, as derived by iGuo et all (|2010h . This SMF ranges 
from stellar masses of 10 s ' 3 Mq to W 117 Mq. In order to ex- 
tend the SMF to even lower stellar masses, we include the 
measurements bv lBaldrv et alj (|2008l ) between 10 7A Mq and 
10 83 M Q . 

As we aim at constraining the stellar mass content of 
haloes back to z — 4 we require high-quality observed SMFs 
at high redshift. In general wide high-redshift surveys cap- 
ture many massive galaxies but are not deep enough to con- 
strain the low-mass tail of the SMF. Conversely, deep sur- 
veys cover too small volumes to constrain the massive tail. 
In practice obtaining a SMF that is equally well constrained 
for massive and low-mass galaxies is difficult, as this would 
involve a survey that is both wide (capturing many mas- 
sive galaxies) and deep (detecting even galaxies of very low 
mass). 

We use the ma s s fun ctions of the wide survey of 
IPerez-Gonzalez et al] (|2008h to constrain the massive end. 
These SMFs have been derived in 12 redshift bins from a 
sample of galaxies at < z < 4 selected in three differ- 
ent fields at 3.6 — 4.5 {im with the Infrared Array Cam- 
era jFazio et al. 2004 ) on board the Spitzer Space Tele- 



scope (|Werner et alj|2004h . The survey covers a total area 



of 664 arcmin 2 and the obtained galaxy catalogue is 90 per 
cent complete down to ~ 5 ^iJy at 3.6 fim and ~ 4 (ijy at 
4.5 /im. Due to the wide area of the survey, the SMFs put 
tight constraints on the abundances of massive galaxies up 
to z = 4, allowing for detailed modeling of the evolution of 
such galaxies. 

For the low-mass e nd we use SMFs in 6 redshift bins 
by ISantini et all (|2012l ). which have been obtained from 
Early Release Science observations taken in the Y, J, 
and i/-bands with the near-IR Wid e Field Camera 3 on 
board HST dWindhorst et all 1201 J ), and from deep K- 
band images taken wi th the near-IR VLT imager Hawk-I 
l|Castellano et al1l201Cll ). The galaxy catalogue has been se- 
lected in the if-band, covers an area of 33 arcmin 2 , and is 
complete to a typical magnitude of K = 25.5. This very 
deep, albeit narrow, survey is able to resolve galaxies of low 
mass up to z — 4.5, and allows us to study the evolution 
of low-mass galaxies in detail. For all SMFs we have con- 
verted the stellar masses to agree with a Chabrier IMF. An 
overview of the 19 SMFs used in this work and their redshifts 
is shown in Figure [6] 
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Figure 1. Upper panel: Sketch of the stellar-to-halo mass ratio 
as a function of halo mass peaking around the characteristic mass 
Mi where it has the normalization N . It has a low-mass slope /3 
and a high mass slope —7. Lower panel: Sketch of the stellar-to- 
halo mass relation as a function of halo mass. The low-mass slope 
is 1 + ft and the high mass slope is 1 — 7. 



iMoster et ahl (|20ld ): 
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+ 



M 
Mi 



(2) 



It has four free parameters: the normalization of the SHM 
ratio N, a characteristic mass Mi, where the ratio is equal 
to the normalization N, and two slopes ft and 7 which in- 
dicate the behavior of m/M at the low and high-mass ends, 
respectively. Figure [T] illustrates the SHM relation in two 
forms, indicating the four parameters. The halo mass where 
the SHM ratio is maximal is related to the characteristic 
mass Mi by 



3 CONNECTING GALAXIES AND HALOES 

In this section we describe how we derive the relationship 
between the stellar mass of a galaxy and the mass of its dark 
matter halo. For this we first need to specify the SHM ra- 
tio. A direct comparison of the halo mass function and the 
stellar mass function shows that the baryons are converted 
into stars with very different efficiencies in haloes of differ- 
ent mass. Star formation is most efficient in systems like 
the Milky Way with M v i r ~ 1O 12 M0, and is substantially 
less efficient in much more massive or much less massive 
haloes. We adopt the parameterization of the SHM ratio of 



M max = Mi 



1/(0+7) 



(3) 



This SHM ratio can then be used to compute the aver- 
age stellar mass in a dark matter halo. However, we expect 
that two haloes of the same mass M may harbor galaxies 
with different stellar masses, since they can have different 
halo concentrations, spin parameters, and merger histories. 
Therefore we assume an intrinsic scatter around this aver- 
age relation when assigning the stellar mass to an individual 
galaxy: we draw a stellar mass m from a lognormal distri- 
bution with a mean value given by equation ([2| and with a 
scatter of <r m (logm) = 0.15. 
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Figure 2. Evolution of the SHM relation parameters with redshift in a model without observational mass errors. The symbols correspond 
to the values that have been derived with the classical abundance matching approach at individual redshifts. Different colors represent 
the different SMFs that have been used to derive the SHM relation: red crosses for the SDSS SMF, green diamonds for the PG08 SMFs 
and blue triangles for the S12 SMFs. The solid line corresponds to a multi-epoch abundance matching model that takes into account 
that satellites are accreted at different epochs. The shaded area indicates the la confidence levels. For Mi and N we assume a second 
order polynomial in z and for /3 and 7 a power law in z. 



3.1 The evolution of the stellar-to-halo mass 
relation 

As a first step, we investigate how the parameters of the 
SHM relation evolve with redshift. For this we assume that 
at a given redshift the relation between the stellar mass of 
a satellite galaxy and the maximum mass of its dark mat- 
ter halo over its history is the same as the SHM relation of 
central galaxies. This assumption is only an approximation 
as the stellar mass of satellites is related to the halo mass 
at infall and the SHM relation is expected to have changed 
since this epoch. However, at any redshift the number of 
subhaloes of a given mass is lower than the number of main 
haloes of the same mass, so the approximation is not greatly 
in error. With this approach we are able study how the pa- 
rameters of the SHM relation depend on redshift, and which 
functional form may be appropriate. 

In order to constrain the parameters at a given red- 
shift (where an observed SMF is available), we populate the 
haloes in the simulation with galaxies using the SHM rela- 
tion to derive their stellar masses. Once the simulation box 
is filled with galaxies, it is straightforward to compute the 
model SMF. We then assign a likelihood to this model SMF 



by relating it to the observed SMF at this redshift: 
C = exp (-xl) 

2 _ 1 ^ / log$ mod (mj) - log$ obB (mi) \ 2 
Xr N *^\ <Tob.(mO ) ■ [> 

Employing a Markov chain Monte Carlo (MCMC) method, 
we sample the probability distribution for the parameters 
and extract the best-fit values and their la errors. We repeat 
this procedure for every observed SMF available and plot the 
results in Figure (symbols with errorbars). 

The transition mass Mi increases with increasing red- 
shift and reaches a maximum around z ~ 2, after which it 
stays constant or even slightly decreases again. This means 
that at z = the conversion of baryons into stars has been 
most efficient in haloes with a virial mass of log(M v i r /Mo) ~ 
11.7, while at higher redshift the conversion was most effi- 
cient in more massive haloes with log(M v ir/A'/Q) ~ 12.3. 
The normalization N, i.e. the SHM ratio at the most effi- 
cient halo mass, decreases with increasing redshift, reaches 
a minimum around z ~ 2 and then increases again. We note 
that this does not imply that individual galaxies lose mass 
as halo mass generally increases for an individual system. 
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Figure 3. Stellar mass of the central galaxy as a function of the 
virial mass of the main halo for different redshifts for the model 
without observational mass errors. The thick lines correspond to 
the ranges where the relation is constrained by data, while the 
thin lines are extrapolations. The arrow shows the evolution of a 
system with Af v j r (z = 0) = 10 14 Mq from z = 1 to 0, indicating 
that if observational mass errors are not accounted for, massive 
systems lose significantly more stellar mass than expected from 
stellar evolution models. 



The low- mass slope ft is only constrained by the S12 
SMFs. This is a consequence of the properties of the surveys 
from which the SMFs are obtained: The survey in PG08 does 
not provide strong constraints on the low-mass slope of the 
SMFs and consequently on /?. The survey used in S12, on 
the other hand, is able to detect low-mass galaxies. From 
these constraints, we see that /3 is decreasing with redshift, 
which means that the slope of the SHM relation becomes 
shallower with increasing redshift. For the high mass slope 
7 constraints are provided only by the PG08 SMFs. These 
constraints indicate that 7 is also decreasing with redshift, 
i.e. the slope of the SHM relation becomes steeper with in- 
creasing redshift. 

3.2 The multi-epoch abundance matching model 

As we now have a first idea how the SHM parameters should 
evolve with redshift, we can construct a self-consistent 
model. As the next step in this multi-epoch abundance 
matching (MEAM) model, we determine a convenient red- 
shift parameterization of the SHM parameters. For the char- 
acteristic mass Mi and the normalization N, we use a second 
order polynomial in z, while for the slopes /3 and 7 we use 
a power law: 



log Mi (z) 
N(z) 

7(2) 



M10 + M11 z + M 12 z 2 
JV10 + N llZ + N 12 z 2 , 
Mz + l)^, 
7io(z + l) 711 . 



(5) 
(6) 
(V 
(8) 



a function of the virial mass of its dark matter halo M v i r 
and redshift z. For satellite galaxies, the stellar mass is now 
determined by employing the SHM relation at the time at 
which it still was a central galaxy, i.e. at the redshift of infall 
Zinf- We then assume that the stellar mass of a satellite 
does not change after its subhalo entered the main halo, 
i.e. we neglect stellar stripping and star formation in the 
satellite. We thus get a redshift- dependent SHM relation 
m(M, z) that depends on the virial mass of the main halo 
and redshift at the current epoch for central galaxies and on 
the infall mass of the subhalo and infall redshift for satellite 
galaxies. The stellar mass of a galaxy at redshift z is then 
given by 



m(M v i r , z) 
m(M inf , z inf ) 



for host haloes 
for subhaloes . 



(9) 



In this way we can derive the stellar mass of any galaxy 
just by employing the redshift-dependent SHM relation for 
central galaxies. 

This relation can then be constrained with a set of SMFs 
at different redshifts. To this end we employ the simulation 
at all redshifts where an observed SMF is available and pop- 
ulate all main haloes and subhaloes in the simulation with 
galaxies using the redshift-dependent SHM relation to de- 
rive their stellar masses. We then compute the N z model 
SMFs for each redshift and assign a likelihood to the model 
by relating them to the corresponding observed SMFs: 



2 



Cxp ^-Xr 



N : 



2 



i=l 



1 y^v" /log«l> mo d(mj) - lag$ oh3 (mj, 



3=1 



(j b s (raj ) 



(10) 



The stellar mass of a central galaxy is then given as 



We again employ an MCMC method, sample the probability 
distribution for all the parameters and extract the best-fit 
values and their la errors. We show the results in Figure [2] 
for the simultaneous fit to all SMFs (black line and shaded 
area). For all four SHM parameters, the lines for the simul- 
taneous fit follow the trends that have been seen for the fits 
at individual redshifts. 

The SHM relation for a set of redshifts between z — 
and z = 4 is plotted in Figure [3] At a fixed halo mass at the 
high-mass end (log(M v i r /MQ) > 13.0), the stellar mass in- 
creases with with increasing redshift. At intermediate halo 
masses (log(Af v i r /M0) ~ 12.0) the stellar mass decreases 
with increasing redshift, reaches a minimum at z ~ 2 and 
then increases again (reflecting the evolution of the normal- 
ization N). The fact that at a given halo mass the mass of 
the central galaxy is lower at present than at an earlier epoch 
does not necessarily imply that individual galaxies lose mass 
during their evolution, as haloes also become more massive 
over time. At a fixed halo mass of log(M v i r /M0) = 12.0 the 
model stellar mass is log(ra*/M0) = 10.5 at z = 4 and only 
log(m»/M0) = 10.2 at z — 2. However, a typical halo with 
log(Mvir/M0) = 12.0 at z — 4 will increase its virial mass 
up to log(M vir /M ) = 12.7 at z = 2, for which the model 
stellar mass is log(ra*/M0) = 11.0, so that the stellar mass 
of the central galaxy does not decrease, but increase. 
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Figure 4. Evolution of the SHM relation parameters with redshift for the model including observational mass errors. The symbols of 
different color correspond to the classical abundance matching approach at individual redshifts: red crosses for the SDSS SMF, green 
diamonds for the PG08 SMFs and blue triangles for the S12 SMFs. The solid line corresponds to the multi-epoch abundance matching 
model. The shaded area indicates the \a confidence levels. For all four parameters we assume a linear dependence on the expansion 
factor a = (1 + z) -1 . 



However, we encounter a severe problem for massive 
galaxies at low redshift: the growth of halo mass cannot 
compensate for the decreasing stellar mass. A typical mas- 
sive halo of log(M V ir/iWo) = 14.0 at z — has grown from 
a z = 1 halo with a virial mass of log(M v ir/Mo) = 13.6. 
The model stellar mass for the central galaxy of this sys- 
tem at z = 1 is log(m*/M0) = 11.4, while at z = it is 
only log(m*/M0) = 11.3. This would imply that massive 
systems lose stellar mass at low redshift and is unphysical, 
as stellar mass loss rates from supernovae are much lower 
than accretion rates onto the central galaxy. This indicates 
that an important effect has been neglected in our model 
and that the effects of this missing ingredient are most se- 
vere at the massive end. In the following section we show 
that neglecting errors in the observed stellar mass impacts 
the results for massive galaxies strongly and that taking this 
into account results in a model where galaxies do not lose 
stellar mass during their evolution. 



3.3 Introducing observational mass errors in the 
model 

Estimating stellar masses for an observed sample of galaxies 
involves obtaining the luminosities of the galaxies and then 
converting these to stellar masses through a stellar popu- 
lation synthesis model and dust corrections. The errors in- 
troduced in all these steps are typically of the order of the 
uncertainties in the stellar populations synthesis analysis. 
Assuming that these errors are Gaussian, the shape of the 
SMF implies that at a given observed mass there are more 
galaxies with a lower true mass than galaxies with a higher 
true mass, i.e. more galaxies are scattered in than scattered 
out. As a result the observed SMF is shallower than the 
true SMF, particularly at the massive end. It is important 
to include this effect in our model. 

We do this by adding a Gaussian scatter to the stellar 
mass of each model galaxy. An 'observed' stellar mass in our 
model is thus obtained in three steps: We first compute the 
average mass for a given halo mass and redshift (or infall 
halo mass and infall redshift for subhaloes), then we add a 
fixed intrinsic scatter of a m (log m) = 0.15 dex to account for 
different halo parameters and histories, and finally we add 
an observational scatter to account for errors in the obser- 
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Table 1. Fitting results for the redshift-dependent stellar-to-halo mass rela- 
tionship. 



Mio 


Mn 


Nio 


Nn 


PlO 




TiO 


711 


Best fit 11.590 
cr 0.236 


1.195 
0.353 


0.0351 
0.0058 


-0.0247 
0.0069 


1.376 
0.153 


-0.826 
0.225 


0.608 
0.059 


0.329 
0.173 



Notes: All masses are in units of Mq . 



vational estimate of stellar mass. The value of this scatter is 
(j obs =0.1 dex for the SDSS SMF (as estimated by LW09), 
(j b s = 0.3 dex for the SMFs by PG08 and S12 at z < 3 and 
(Jobs = 0.45 dex for the SMFs at z ^ 3 (see appendix B of 
PG08). 

3.3.1 The evolution of the SHM parameters 

First we employ the classical abundance matching approach 
as described in Section |3~T1 now taking into account observa- 
tional mass errors, to investigate how the SHM parameters 
evolve with redshift, and if a different functional form has 
to be adopted. The results are plotted in Figure 0] (sym- 
bols with errorbars). As for the fit without observational 
mass errors, the transition mass Mi increases with increas- 
ing redshift. This increase is strong at low redshift, while 
at high redshift Mi increases only slightly. Unlike before 
we do not find that Mi has a maximum beyond which it 
decreases again. This result implies that at high redshift, 
massive haloes with log(M v i r /M0) ~ 12.5 are most efficient 
in converting baryons into stars, while at lower redshift the 
halo mass for which the baryon conversion efficiency is max- 
imal, decreases towards log(M v j r /M0) ~ 11.6. 

The normalization N (the conversion efficiency at the 
peak) is now found to decrease monotonically with redshift. 
The reason why the normalization does not change much 
at low redshift with respect to the model without observa- 
tional scatter, but is very different at high redshift, is the 
evolution of the low-mass slope a of the SMF. At low red- 
shift this slope is very shallow, so relatively few galaxies are 
scattered from lower masses towards the characteristic mass 
of the SMF (which corresponds to the peak efficiency). At 
high redshift, the low-mass slope is much steeper, so that 
more galaxies are scattered towards (and beyond) the mass 
of the peak efficiency. As a result, the observed SMFs at low 
redshift are very similar to the true SMFs, while at high red- 
shift the observed SMFs are shifted towards higher masses. 
Therefore it appears as if the normalization is increasing 
at high redshift, if the observational mass error is not taken 
into account. When we account for this, we find that the true 
SMFs are shifted to lower masses and the normalization of 
the SHM ratio N no longer increases at high redshift. 

The results for the low-mass slope /3 are very similar 
to what we found before. This is a consequence of the shal- 
low low-mass slope of the SMFs, so that the observational 
mass error does not change the SMFs much. As before, /3 is 
decreasing with increasing redshift, so that the slope of the 
SHM relation becomes shallower with increasing redshift. 
Again, f3 is only constrained by the S12 SMFs. 

The high-mass slope 7 is very different to what was 
found without observational mass errors. Where before 7 



was decreasing with increasing redshift, we now find that 7 
is increasing. The reason for this is the steep slope of the 
SMFs at the massive end, which means that the true SMFs 
are even steeper than the observed SMFs. Therefore the high 
mass slope of the SHM relation has to become shallower 
which yields a larger 7. Here, the constraints are provided 
only by the PG08 SMFs, as only this survey is wide enough 
to detect very massive galaxies. 

3.3.2 The multi-epoch abundance matching model 
including observational mass errors 

Now that we have checked how the SHM parameters evolve 
with redshift, we can again employ our self-consistent model. 
As this evolution is found to be different with respect to the 
model without observational mass errors, we first need to 
determine the new redshift parameterization of the SHM 
parameters. As all parameters evolve considerably at low 
redshift, but evolve only slightly at high redshift, we found it 
more convenient to use a linear dependence on the expansion 
parameter a — (z + 1) : 



log Mi (2) 


= Mio + Mn(l 


-a) =Mio + Mn— ^-p 

2+1 


,(H) 


N(z) 


= JVio + JVn(l- 


-0) = Nio+Nu-^—, 
z + l 


(12) 


P(z) 


= fto + /3ii(l- 


a)=p w +Pu , 

Z + 1 


(13) 


7(*0 


= 710+711(1 — 


. Z 

a.) = 710 + 711 — rr- 
z + l 


(14) 



The stellar mass is given by equation ((2| as a function 
of the virial mass of the main halo Af v i r and redshift z for 
central galaxies and as a function of the infall mass of the 
subhalo Minf and infall redshift Zi n f for satellite galaxies (see 
equation ©). We constrain the eight parameters Mio, Mn, 
Nio, A^n, Pio, fin, 710 and 711 by fitting our model to a set of 
SMFs at different redshifts. Again we employ the simulation 
at those redshifts where an observed SMF is available and 
populate all main haloes and subhaloes in the simulation 
with galaxies using the redshift-dependent SHM relation to 
derive their stellar masses. We then compute the 19 model 
SMFs and employ an MCMC method to find the best-fit 
values and their la errors. The results are summarized in 
Tabled] 

We show the results in Figure [4] for the simultaneous 
fit to all SMFs (black line and shaded area). The redshift 
evolution of all four SHM parameters in the simultaneous 
fit follow the trends that have been seen for the fits at indi- 
vidual redshifts. We plot the SHM relation and the baryon 
conversion efficiency for a set of redshifts between 2 = and 
z = 4 in Figure [5] At a fixed halo mass, the stellar mass now 
decreases with increasing redshift for low mass, intermediate 
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Figure 5. Stellar-to-Halo mass relation for central galaxies as a function of redshift. This relation is also valid for satellite galaxies when 
virial mass and current redshift are substituted with subhalo infall mass and infall redshift, so that the satellite stellar mass is obtained 
when it was last a central galaxy and does not change until it is accreted onto the central galaxy. Left panel: Evolution of the SHM 
relation. Right panel: Evolution of the baryon conversion efficiency m* iCGn /Mj, = fn,*,con/(Af v ; r /i,) = (mt, iCen /M v j r ) X (f2 m /f2(,). The 
thick lines correspond to the ranges where the relation is constrained by data, while the thin lines are extrapolations. 



and massive galaxies. As a consequence, all central galaxies 
grow with time. Moreover, independent of halo mass, we find 
that in all systems the central galaxies grow faster in stellar 
mass than the respective main haloes accrete dark matter. 

Finally we compare our model SMFs to the observed 
ones in Figure [U The model SMFs have been computed 
using the result of the simultaneous fit to all SMFs, i.e. 
the parameters given in Table [1] The dashed lines show the 
'true' SMFs, that have been computed after intrinsic scatter 
has been added to the galaxy catalogue, but before obser- 
vational mass error. The solid lines show the model SMFs, 
after mass errors have been included. At low redshift and at 
low stellar masses, the observational mass errors do not im- 
pact the shape of the SMFs significantly, as the slope of the 
SMF is shallower than at higher redshift and at the massive 
end. Beyond the characteristic stellar mass, the 'true' SMFs 
are much steeper than the observed ones, so that there are 
substantially fewer very massive galaxies at high redshift 
than the observational SMFs appear to indicate. 

The model SMFs including observational mass error are 
in very good agreement with the observed SMFs. This in- 
dicates, that our parameterization of the redshift evolution 
of the SHM relation is compatible with observational con- 
straints. However, there are a number of discrepancies, that 
stem mainly from the fact that the observed SMFs of dif- 
ferent authors are not in agreement. The SMF of PG08 at 
low redshift (z — 0.1) is higher than our model, while the 
SDSS SMF is fit very well. The reason for this is that the 
SDSS has much smaller errors than the PG08 SMF, so that 
during the fit a much higher likelihood is given to the model 
that agrees with the SDSS. At higher redshift there are dis- 
crepancies between the model and the observed SMFs at the 
low-mass end. While the model overpredicts the abundance 
of low-mass galaxies with respect to the PG08 SMFs, it un- 
derpredicts the abundance with respect to the S12 SMFs. 
This is due to the PG08 SMFs having a much shallower 



slope than the S12 SMFs. As the errors on both sets of ob- 
servations are of the same order, a model is favored that 
yields average abundances. 



4 THE EVOLUTION OF THE STELLAR 
CONTENT 

Having constrained the evolution of the SHM relation for 
central galaxies, we can now use our model to infer the evolu- 
tion of the stellar content of galaxies. To this end we employ 
the dark matter halo merger trees that have been extracted 
from the MS and MS-II. These merger trees connect haloes 
at any given epoch to their descendants at later epochs 
or to their progenitors at earlier epochs. As our redshift- 
dependent SHM relation determines the connection between 
galaxies and haloes at any epoch, we can use the connection 
between haloes across time to infer the connection between 
galaxies across time. By averaging over a sample of haloes 
with the same mass at z — 0, our model yields the average 
stellar mass growth of central galaxies as a function of halo 
mass. We apply this method for sets of halo masses between 
1O 11 M0 and 10 14 Mq. This analysis is rela ted to previous 
studies that employed semi-ana l ytic models ( Dc Lucia et al.l 
2006); [5e Lucia fc Blaizotj[2007l ; lGuo fc White! l20Qg| ). Fora 
comparison and discussion we refer to Section [5] 

As a first step, we select a sample of merger trees by 
fixing the virial mass of the main halo at z = 0. In the next 
step, we populate the haloes in these merger trees using 
our redshift-dependent SHM relation to compute the stellar 
masses of their central galaxies. At every epoch, i.e. simu- 
lation snapshot, we then compute the average stellar mass 
of the central galaxies in the main haloes. Finally, we calcu- 
late the difference of the average central galaxy stellar mass 
between two subsequent snapshots and divide this by the 
time that has elapsed between those snapshots to derive the 
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Figure 6. Comparison of the model (lines) and the observed SMFs (symbols). The redshift is indicated in each panel. The model SMFs 
have been computed using the results of the fit including observational mass error. The solid lines show the SMFs when mass errors are 
included and represent a fit to the observations. The dashed lines show the 'true' SMFs, before observational mass error is added. 



average stellar mass growth rate of the central galaxies. By 
doing this for a set of trees in halo mass bins, we can com- 
pute the average growth rate as a function of the z — halo 
mass and redshift. 

In the left panels of Figure [7] we show the average frac- 
tion of 2 = mass assembled as a function of redshift for 
dark matter haloes and central galaxies. For four 2 = 
halo masses, the assembly history of a central galaxy is 
compared to that of its parent halo. In low-mass systems 
with log(M v i r /MQ) < 12.5 the halo grows much slower than 
its central galaxy. This indicates that in these systems the 
growth is dominated by internal gas physical processes such 
as cooling and star formation. Only rather little external 
material is falling into the halo. For massive systems with 
log(A4vir/Me) > 12.5 the trends are opposite at low red- 
shift: the fractional growth of haloes since z ~ 2 is much 
stronger than the stellar growth of their central galaxy. A 



large fraction of the stellar mass of massive galaxies is al- 
ready in place at z < 2. In these massive systems the growth 
is dominated by external processes at late times, as they are 
accreting large amounts of material. While the halo grows 
rapidly, the central galaxy does not, as most of the accreted 
systems take a very long time to reach the centre of the halo 
and instead survive to z — as 'satellite' members of the 
group or cluster. 

To further study when most of the stellar and halo 
mass was in place we fol low and extend the analysis by 
IConrov fc Wechslerl (20091 ) . We plot the average assembly 
redshift of main haloes and central galaxies as a function of 
2 = mass in the right panel of Figure [7| We show the red- 
shift at which 25, 50, and 75 per cent of the present stellar 
and halo mass was already in place. In low-mass systems, 
the stellar mass assembled much later than the halo mass. 
For a system with M v i r = 1O 12 M0, half of the stellar mass 
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Figure 7. Left panels: Average fraction of z = mass assembled as a function of redshift for dark matter haloes (blue dashed lines) 
and central galaxies (red solid lines). Each panel compares the mass assembly history of a central galaxy to that of its parent halo. The 
z = halo and stellar masses are indicated in each panel. While for low-mass dark matter haloes most of the mass assembles at early 
times, massive haloes only assemble late. For galaxies these trends are opposite. Right panel: Average formation redshift of dark matter 
haloes (blue dashed lines) and central galaxies (red solid lines) as a function of z = mass. The three different lines indicate the redshift 
at which 25, 50, and 75 per cent of the mass was in place. 



of the central galaxy forms after redshift z = 0.7 while half 
of the virial mass was already assembled by redshift z = 1.3. 
In more massive systems haloes assemble later while galax- 
ies assemble earlier: for a system with a z = halo mass 
of M v i r = 10 13 Mq, half of the halo has assembled only by 
redshift z — 1.0, but half of the stellar mass was already in 
place by redshift z — 1.8. 

Very massive galaxies in haloes with Af v i r = 1O 14 M0 
show a very interesting behavior: while they begin to form 
very early and have 25 per cent of their stellar mass in place 
already by redshift z = 3.0, it takes them a long time to fully 
assemble. Half of their stellar mass is in place only by red- 
shift z — 1.6, while 75 per cent is assembled by z = 0.3. The 
final z = stellar mass is then assembled very quickly. The 
reason why these massive galaxies grow very fast at early 
times, only grow slowly at intermediate redshifts and grow 
fast again at late times can be found in the two processes 
that contribute to the growth of galaxies. At high redshift, 
a massive galaxy grows effectively by star formation, while 
at low redshift, stellar accretion leads to a fast growth. At 
intermediate redshift however, neither of these processes is 
effective, so that the galaxy grows only very slowly. In the 
following section we aim to constrain and separate these two 
processes. Note that while massive galaxies assemble their 
final stellar mass only late, the stars in these galaxies have 
formed much earlier, such that this result is not in conflict 
with the observation that the stars in massive systems are 
older than stars in low-mass galaxies. 



4.1 Star formation and accretion histories 

The stellar mass of a galaxy can grow either by star forma- 
tion inside the galaxy (in-situ) or by the assembly of stars 



that formed outside the galaxy and have been accreted later 
(ex-situ). In principle, both processes can contribute simul- 
taneously to the growth of a central galaxy, but in practice 
one process will dominate over the other as a function of halo 
mass and redshift. In order to separate the contributions 
from star formation and accretion, we have first to deter- 
mine the amount of accretion of subhaloes onto the central 
object in the halo. Using our model, we can then populate 
these subhaloes with satellite galaxies and compute the stel- 
lar accretion rate. Finally, the SFR can be determined from 
the difference between the total growth rate and the accre- 
tion rate. 

We use the same set of dark matter halo merger trees 
as before, i.e. we use all trees with a common z = virial 
mass of the main halo and assign a stellar mass to all haloes 
and subhaloes using our redshift- dependent SHM relation. 
For every simulation snapshot, we then record all satellite 
galaxies that are accreted onto the central galaxy until the 
time of the next snapshot. We note, that, by construction, 
this also includes orphan galaxies. At each snapshot we then 
sum the stellar mass of all recorded merging galaxies and 
divide this by the number of central galaxies in the sam- 
ple (which corresponds to the number of merger trees). We 
further assume that a fraction of satellite stars / csc = 0.2 
escapes from the central galaxy and ends up in the halo 
as diffuse stellar material not detected in standard surveys. 
This fiducial value is motivated by hi gh-resolution simu- 
lations of binary galaxy mergers (e.g. Nipoti et alj 120031 : 
iNaab fc Burkertll2003l : iMoster et at] l201ll ). For a different 
value the resulting stellar accretion and star formation rates 
change accordingly (e.g. for / osc = the accretion rates are 
25 per cent higher). In this way we get the average stellar 



12 B. P. Moster et al. 



^ 10 










in — situ 






ex — situ 

total 
















50% 








10% 











M 
O 



10 

9 

8 
10° 

10" 1 

10" 8 

10" 3 



1 2 3 4 6 
z 



2 3 4 6 
z 



1 2 3 4 6 
z 



1 2 3 4 6 
z 



Figure 8. Star formation and accretion histories for different z = halo masses. The rows from top to bottom show the star formation 
history, the stellar accretion history, the stellar mass formed in-situ and ex-situ, and the fraction of ex-situ formed stars. Each column 
represents a different z = halo mass. In the first two rows the black lines and shaded areas represent the rates that have been derived 
with our model and their lcr confidence levels. The red lines are fits to these relation as given by equations lll7t - H20H and H21IH23I I. The 
blue line in the first row shows the stellar mass loss rate. In the third row, the red dotted lines give the stellar mass of in-situ formed 
stars that have survived, the blue dashed lines give the stellar mass that has formed ex-situ, and the black line is the total stellar mass 
in the central galaxy. The dotted lines in the fourth row indicate that 50, 10 and 1 per cent of the stellar mass has formed ex-situ. 



accretion rate onto the central galaxy as a function of the 
z = halo mass and redshift. 

Having determined the total stellar growth rate and the 
amount of growth that is due to accretion, the SFR of an av- 
erage galaxy can be computed from the residual stellar mass 
growth attributed to star formation. However, the relation 
between this residual growth rate and the SFR is compli- 
cated by mass loss due to dying stars: 

dm rcs _ dm to t dm acc _ dmsF dmi oss _ dmis , . 
dt ~ ~dl df~ _ dt df~ _ dt ' 

This equation states that the residual growth rate dm r es/dt, 



which is defined as the difference between the total growth 
rate dm to t/dt and the accretion rate dm aC c/dt, is given by 
the difference between the SFR dmsF/di and the stellar 
mass loss rate drai oss /di, and corresponds to a net in-situ 
growth rate dmis/dt. It can be solved iteratively, starting 
with a SFR that is equal to the residual growth rate at every 
epoch. We go to the time of the first snapshot and compute 
the amount of stellar mass that is formed in an average 
galaxy. Then we move to the second snapshot and compute 
both the stellar mass that forms at this time and the stellar 
mass that is lost in the population that formed in the first 
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snapshot. For the latter we employ the formula 

/-W = 0.Q5xln(i±f^l) , (16) 

where /i oss (t) is the fraction of mass lost by a stellar popula- 
tion with an age t. This formula is a fit t o the mass loss of a 
ste llar population with alChabrierl (|2003l ) IMF as developed 
bv lBruzual fc Charlotl (|2003T ) and implies that 40 per cent 
of stellar mass is lost within 1 Gyr and 50 per cent of stellar 
mass is lost within 10 Gyr. 

We then move to the third snapshot, compute the stel- 
lar mass that forms at this time and derive the stellar mass 
that is lost in the populations that formed in the first and 
second snapshots. This procedure is continued until z = 0, 
while at every snapshot the amount of lost stellar mass is 
the sum of the lost mass of all populations that have formed 
until then. We can then compute the loss rate between subse- 
quent snapshots and the in-situ growth rate as the difference 
between SFR and loss rate. This in-situ growth rate is then 
compared to the residual growth rate as described by equa- 
tion [15] The difference between those rates is then added to 
the SFR (separately at each snapshot) and the procedure is 
repeated and iterated until the in-situ growth rate and the 
residual growth rate are equal. We repeat this procedure for 
each sample of halo merger trees with a given z — halo 
mass. 

The resulting SFRs and accretion rates are shown in 
the first and second rows of Figure [8] (black solid lines), 
respectively, for four different z = halo masses. The la 
confidence levels for both rates (shaded areas) have been 
computed by employing the upper and lower estimates for 
mean stellar masses according to our SHM relation. We have 
also added the stellar mass loss rates in upper panels (blue 
line). All results represent the average rates, i.e. the average 
of both red and blue galaxies at a given z — halo mass. 
Of course, individual galaxies can have higher or lower SFRs 
depending on other galaxy properties like color, so that blue 
(red) galaxies can have higher (lower) SFRs. 

We find that the star formation history of central galax- 
ies peaks at a certain redshift. This peak SFR redshift de- 
pends on halo mass: for massive systems this peak occurs 
at high redshift (z = 3 for M vil = 1O 14 M ) while for 
lower masses the SFR peaks at a lower redshift [z = 0.5 
for M v i r = 1O 12 M0). For systems with an even lower mass 
of M v i r J$ 1O 11,7 M0 the peak has not been reached yet by 
2 = 0. After the peak has been reached, the SFR quickly 
drops, such that massive galaxies are quenched at z ~ 2. 
Furthermore, the peak SFR is higher in more massive sys- 
tems. We also find that in systems that are more massive 
than M v i r = 10 12 ' 5 Mq, the stellar mass loss rate due to dy- 
ing stars is larger than the SFR at z — 0, such that the mass 
of stars present in the central galaxy is decreasing. 

The accretion rates increase with time for all systems. 
Massive systems have a high accretion rate that dominates 
the total growth rate at late times, while for low-mass sys- 
tems, accretion is insignificant at all times. For a system 
with Mvir = 1O 14 M0 the accretion rate dominates over the 
SFR at 2 < 1. The halo mass at which the accretion rate and 
the SFR are equal at z = is M v i r = lO 12 ' 7 A/ . For lower 
masses the contribution of accretion to the total growth rate 
quickly drops. The reason why accretion contributes signif- 
icantly to the overall growth in massive systems and it is 



insignificant in low-mass systems can be inferred from the 
SHM relation. Massive haloes host relatively massive sub- 
haloes, and these subhaloes are filled with relatively massive 
satellite galaxies. The subhaloes that have e.g. one tenth of 
the main halo mass host galaxies that are almost as massive 
as the central galaxy, so that the resulting accretion rates 
are high. In low-mass systems, however, the subhaloes are 
almost entirely devoid of stars, so that the accretion of the 
subhaloes onto the center does not result in a significant 
stellar accretion rate. Moreover, Figure [7] shows that low- 
mass haloes have largely assembled by z — 1, so that any 
late time growth has to come from star formation. 

To compare further the relative importance of star for- 
mation and accretion for systems of different mass, we plot in 
the third row of Figure[8]the amount of stellar mass that has 
formed in-situ, i.e. through star formation, and the amount 
that has formed ex-situ, i.e. has been accreted. In the fourth 
row we show the fraction of stellar mass of the central 
galaxy that has formed ex-situ. In massive systems with 
Mvir ^ 1O 13 M0, all of the stars that have formed in-situ and 
survived were already in place at high redshift (z > 2). Al- 
though, the loss rate is higher than the accretion rate at late 
times, and the amount of in-situ stellar mass has decreased 
since z ~ 2, the total stellar mass of the system increases 
with time, as the accretion rate compensates for the lost 
stellar mass. In very massive systems with M v i r = 10 14 Mq 
more than half of the central stellar mass is formed ex-situ. 
In lower mass systems the amount of ex-situ stellar mass is 
decreasing considerably: for M v i r = 10 13 Mq the fraction of 
ex-situ stars is 20 per cent, for M v i r = 10 12 Mq the frac- 
tion is 2 per cent, and for M v i r = 10 12 M© the fraction is 
1 per cent. At high redshift (2 > 2), the amount of stellar 
mass that has formed ex-situ is negligible for systems of all 
masses. We further note that the accretion rates grow more 
rapidly in massive systems than in systems of low mass. 

The connection between star formation and accretion 
history and the virial mass of the system can be of general 
use in constraining models of cooling, star formation and 
feedback in galaxies. These relations have been derived en- 
tirely from observations with an empirical model. Although 
this model makes no reference to the physical processes that 
drive the evolution of galaxies, it can be very helpful in order 
to constrain more physical models. In order to easily com- 
pare our model predictions to the results of other models 
such as simulations and semi-analytic models, we provide 
simple fitting functions to the star formation and accretion 
histories. 

We find that at a given 2 = main halo mass the SFR 
is well fit by the relation: 

= / 1 aA«p(A ( i. _a)) 

- ^ 1+ ^ /2exp (fr^)- <"> 

The parameters fi, f%, and fa depend on the 2 = halo 
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mass of the system and follow the relations: 

(log(M vir /M )-/ n ) 2 
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(18) 
(19) 
(20) 



The parameter /1 can be interpreted as the SFR at z — 0, 
/a is the slope of the star formation history at high redshift, 
and h corresponds to the expansion factor a of the peak 
SFR. The values of the fitting parameters are: 



h 



■■ 0.424, 



ho = 2.658 M Q yr" 1 , /n = 11.937, 
ho = 5.507, /21 = 2.437, 

/ 30 = -0.915, /31 = 0.696, / 32 = -0.159 . 

The SFRs obtained with this fitting function are compared 
to the inferred ones in Figure[8](red lines). Overall the global 
fit is very close to the model predictions. For very massive 
galaxies with M v i r = 1O 14 M0 at low redshift (z < 1), the 
fitting function is lower than the model SFR, although it is 
still within the la confidence levels. 

The accretion history for systems of a given z = main 
halo mass can be described with an exponential dependence 
on redshift: 



dm a 



dt 



gi exp 

32 



(21) 



The parameters g\, and gi are functions of the z — halo 
mass of the system and are given by: 
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(22) 
(23) 



The parameter g\ gives the stellar accretion rate at z = 0, 
and <?2 is the slope of the accretion history. The values of 
the fitting parameters are: 



g\o — 1.633 Mq yr 
320 = 0.855, 



3ii 
321 



-2.017, 
0.098, 



312 

322 



-0.806, 
-0.797 . 



The accretion rates obtained with this fitting function are 
compared to the inferred ones in Figure [8] (red lines). The 
global fit is very close to the model prdictions for all redshifts 
and halo masses. 



4.2 Cosmic star formation rate density 

In order to check whether our model SFRs are in agree- 
ment with observational constraints, we compare our model 
prediction for the evolution of the cosmic SFR density 
to data in the literature. For this we compute the SFR 
for each halo in the simulations using our recipe as given 
by equations (|17[l - (|20[ l and compute the volume density 
for each snapshot. We compare the result to observa- 
tional estimates derived with different methods in Figure 
[9] This comparison includes estimates derived from the UV 



9J_rms comparison includes estimates derive d trom tne U V 
llSalim et al ]l2007l; Ivan der Burg et ai1l2010l; Bou wens et al.l 
20091: iRobotham fc Driver! l201ll; iBouwens et all l201ll : 



Cucciati et al.ll2012T i, Ha l|Lv et alj|201ll ). both UV and IR 
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Figure 9. Star formation rate density as a function of redshift. 
The symbols of different colors give observational estimates de- 
rived with different methods. The black line is the result of our 
model. The model was constrained with observed SMFs at dif- 
ferent redshifts and has not been fit to observed SFRs. Still the 
model is in excellent agreement with the data. 



jZheng et al.ll2007i ), the FIR dRuiopakam et al] |201Ch and 
from radio observations (|Smolcic et alj 120091 : iDunne et al.l 
I2OO9I : iKarim et alJl201lh . 

Our model is in excellent agreement with the data at all 
redshifts. We note that our model was only constrained with 
observed SMFs at different redshifts and has not been fit 
to observed SFRs. This means that the observed SMFs are 
broadly consistent with observational inferences for SFRs 
once observational mass errors are taken into account. More- 
over, this shows that our choice for the parameterization for 
the redshift-dependent SHM relation is consistent with ob- 
servational constraints, and that our inferred recipe for the 
mass-dependent star formation histories gives reliable esti- 
mates. 



4.3 Specific star formation rates 

As a further test whether our model predictions are in agree- 
ment within the observational limits, we compute the spe- 
cific star formation rate (SSFR) for galaxies of a fixed stellar 
mass. Observationally, it is very difficult to infer the SSFR as 
a function of the z — stellar mass since observations cannot 
connect galaxies at different epochs. Therefore it is not pos- 
sible to follow the star formation history of an average galaxy 
across time with observations up to high redshift. However, 
it is possible to derive the average SSFR for galaxies within 
a stellar mass range at different redshifts. We can compare 
these observational data to the predictions of our model, by 
populating all haloes in a given snapshot with galaxies us- 
ing our SHM relation to get stellar masses and our recipe 
for the star formation history to get their SFR. We show 
this comparison for low-mass galaxies with log(m»/M v i r ) = 
9.5 in Figure [lOl and for massive galaxies log(ro*/M v ir) = 
10.5 in Figure 1111 The comparison includ e s estimates 
deriv ed from t he UV (jFeulner et al . 2005; Salim ~et al.l 
| 2007|). the IR (lLabbe et all I2OI0I : ISchaer er fc de Barrosi 



2010; Mc Lure et al. 120111 ) , both UV and IR Kheng et al 



20071 : iNoeske et all 120071: iKaiisawa et al] 120101), and from 
radio observation s 1 Daddi et al.l 20071 : Dunne et al. 1 120091 : 
IKarim et al1l201ll ). 
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Figure 10. Specific star formation rate for galaxies with a stellar 
mass of log(m* /M v ; r ) = 9.5 at the various epochs. The symbols of 
different colors give observational estimates derived with different 
methods and the black line is the result of our model. Although 
this relation has not been fit, the model is in excellent agreement 
with the data. 




Figure 11. Same as Figure [lOl but for galaxies with a stellar 
mass of log(m*/M v ; r ) = 10.5 at the various epochs. 



For both stellar mass ranges, our model is in excellent 
agreement with the data, although this relation has not been 
fit. This shows that mass-dependent quantities, such as the 
SSFR can be predicted up to high redshift in addition to 
global quantities such as the cosmic SFR density. As our 
model is able to fit a wide range of observed data, we are 
confident that it can be used to estimate reliably relations 
that cannot be easily observed, such as the star formation 
history of an average galaxy as a function its z — halo 
mass. 



5 DISCUSSION 

In this study we draw conclusions on the assembly histories 
of galaxies over a wide range in mass from a conceptually 
simple approach matching observed galaxy mass functions 
and halo accretion histories from modern high resolution 
dark matter simulations which represent the underlying cos- 
mological model. The early domination of in-situ star forma- 
tion and the growing importance of late accretion of ex-situ 
stars for galaxies residing in more massive haloes found here 
is in good agreement with predictions from semi-analytical 



models (|Kauffmannlll996l ; iKhochfar fc Silkll2006r). Further- 
more, our analysis confirms the results of De Lucia et al.l 
(2006) who find using semi-analytic models that more mas- 
sive (elliptical) galaxies have star formation histories that 
peak at higher redshifts than lower mass systems. 

Although star formation ceases at low redshift for mas- 
sive systems, the stellar mass of their central galaxies can 
grow considerably through accreted satellite galaxies. The 
majority of these ex-situ stars have formed when their galax- 
ies still were centrals, i.e. before they entered the main 
halo. Therefore massive galaxies consist of old stellar pop- 
ulations. On the other hand, massive galaxies start to as- 
semble very early due to star formation, but their final mass 
is a ssembled only very l ate. E mploying semi-analytic mod- 
els, iDe Lucia fc Blaizotl (|2007h find that half of the mass of 
central galaxies of very massive haloes with M V \ T = 1O 15 M0 
forms after z ~ 4, but is locked up in the central galaxy only 
after z — 0.5. This is in good agreement with our results: 
extrapolating the 50 per cent assembly redshift in Figure [7] 
to M vir = 10 15 Af Q roughly yields z = 0.5. 

Our model indicates that low-mass (massive) galaxies 
form late (early), although their dark matter haloes assem- 
ble early (late). Although the haloes of low-mass galaxies 
assembles early it takes a considerable amount of time for 
the gas to cool and fall to the centre and form stars. In 
massive systems, on the other hand, the central galaxies 
(and the galaxies that will become satellites in these sys- 
tems) can start to grow very early and form the bulk of 
their stars at high redshift. However, their haloes assemble 
only late, and the satellites are accreted onto the central 
galaxies at low redshift. Therefore, the stellar populations 
of massive galaxies are old, although massive haloes assem- 
ble late. This demonstrates that 'downsizing' is not anti- 
hierarchi cal, but a consequence of the bottom-up ACDM 
scenario (jNeistein et al.ll2006l ; iFontanot et alj|2009h . 

The contribution of star formation and accretion to the 
total stellar growth of central galaxies depends both on halo 
mass and redshift. While at high redshift accretion is negli- 
gible in all systems, at low redshift star formation dominates 
only in low-mass systems. The halo mass at which both pro- 
cesses contribute equally to the total growth of the central 
galaxy at z — is log(M v i r /Mq) ~ 12.7, which corresponds 
to a stellar mass of log(m/M0) ~ 10.9. Here, both modes 
add ~ IMq yr -1 to the stellar mass of the galaxy. This is 
in good agreement with the results of iGuo fc White! (|2008l ) 
who employ a semi-analytic model to study star formation 
and stellar accretion rates as a function of stellar mass and 
redshift. At z = they find that the stellar mass of cen- 
tral galaxies for which both processes contribute equally is 
roughly log (m/M©) ~ 10.6, with rates of ~ 0.5 -1M© yr -1 
(see their Figure 3). For systems of lower (higher) mass, 
both methods indicate a lower (higher) contribution of ex- 
situ stars. At z = 2 the semi-analytic model predicts that 
the stellar mass of galaxies with an equal contribution of 
the two processes is roughly log(m/Mo) ~ 11.0, with rates 
of - 20M© yr" 1 , while in our model only the most massive 
galaxies with log(m/Mo) > 11.4 have accretion rates that 
are comparable to their SFRs. 

In the present study on average 50 per cent of the stars 
in galaxies in massive haloes (1O 14 M ) at z = are made 
in external galaxies and accreted at later epochs. However, 
in lower mass haloes, like that of the Milky Way (10 



112 Mq) 
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Figure 12. Comparison between central galaxy formation efficiencies found in numerical simulations at different redshifts. Each panel 
corresponds to the indicated redshift. The solid black lines give the average conversion efficiencies needed to fit the observed SMFs and 
the shaded areas indicate the lcr confidence levels . The s ymbols show the result s of hydrodynamical zoom-in simulations run with the 
Gasoline code (blue aster isks: iBrooks et al. 201f, pluses: fGovernato et al .1120121. crosses: Stinson et al. in prep.) and the Gadget code 
(red dots: lOser et al.ll2010l . diamonds: IScannaDieco et al.ll201ll . squares: |Genel et al 1l2012l . triangles: Ipkamotol Ejlj ) . The colored lines 
show the conversion efficiencies predicted by the semi-analytic model bv lGuo et al.N201ll . red dotted line) . While many simulations agree 
well with the predicted conversion efficiency at z = 0, most have too high values at earlier epochs, indicating that they form their stars 
too early. 



this fraction is only a few per cent, making stellar accre- 
tion relatively unimportant for the evolution of the central 
galaxies. This in accord with sta ndard formation and enrich- 
ment scenarios for disk galaxies (iMatteucci fc Francoisll 19891: 



Mo et alj 


1998; 


Prantzos &: Silkl 


1998: iNaab & Ostrikerl 


2006; Kauffmann ct al. 2012). 



Similar trends are reported from three-dimensional cos- 
mological hydrodynamical zoom-in simulations, however, 
with larger fractions of accreted stars for galaxies of all 
masses. For the most massive galaxies this can be up to 
80 per cent (|Oser et al.l 120101 : IFeldmann et al.ll2010l ). Still, 
it is likely that this fraction is reduced to 50 percent (in 
the same mass bin) by the inclusion of strong stellar feed- 
back which reduces the s tellar content of the (typ ically lower 
mass) accreted galaxies (Hirschmann et al. 2012). In zoom- 
in simulations individual stellar particles can directly be fol- 
lowed and it has been shown that e x-situ stars predom i- 
nantly assemble at larger radii (e.g. iLaporte et al.1 12OI2T ). 
In massive spheroidal galaxies, where the fraction of ex- 
situ stars is high, this material plays a significant role in 



the structural evolution of the galaxies (IN aab et al 



Hopkins et aLll2009l ; IFeldmann et al1l2010l : lOser et al 



2009 



2010 



20121 ). In lower-mass disk-like galaxies the accreted stars 



contribute to the build-up of the outer stellar halo and can 
also directly contribute on the order of 10 percent to the stel- 



lar mass of the disk (e.g. lAbadi et al.ll2006l : iFont et al.ll201ll : 
iTissera et "all 1 20121 ). Interestingly, these simulations predict 
ex-situ fractions of more than ~ 20 per cent (up to 40-50 
per cent) which is an order of magnitude higher then found 
in our study. 

Our model constrains the mass budget of cosmological 
galaxy formation simulations not only at the present day 
but also at higher redshift. We compare our model to the 
results of hydrodynamical zoom-in simulations and semi- 
analytic models at six redshifts in Figure 1121 While some 
older simulations have too high conversion efficiencies at z — 
0, most of the newer simulations show good agreement with 
the efficiencies needed to reproduce the local observed SMFs. 
However, almost all simulations have far too high efficiencies 
at earlier epochs. This indicates that these simulations form 
too many stars early (2 ~ 3), while at late times (z < 2) 
star formation is quenched too strongly, so that most of the 
stellar mass at z = is already present at high redshift. 
The lack of late star formation in simulations could possibly 
explain the difficulties in getting prominent discs. Most stars 
are formed at times where processes leading to the formation 
of spheroids are common, e.g. m ajor mergers or spin flips 
|Sales et al.ll201ll : lOkamotol [20121) . If the majority of stars 
form after these epochs, a more prominent disc forms. 

The results of the semi-analytic models agree very well 
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with our relation at z = 0. This is not surprising as these 
models were tuned to fit the local SMF. At higher redshift 
however, their predictions differ from our model, indicating 
as noted in the original papers that they do not yield good 
fits to the observed high redshift SMFs. While the normal- 
ization and the slopes are in good agreement with our rela- 
tion, the halo mass where the efficiency peaks (Mi), moves 
to lower masses with increasing redshift. 

We note that a direct comparison between stellar 
masses measured in simulations and observed stellar masses 
is difficult. Especially at high redshift, it is complicated to 
estimate stellar masses observationally, as many different 
effects have to be taken into account. Moreover, one has to 
make a number of assumptions when fitting the spectral en- 
ergy distributions (SEDs). An alternative way to compare 
stellar masses in models and observations is by making mock 
images of simulated galaxies using a stellar populations syn- 
thesis code and then obtaining a 'photometric' stellar mass 
with an SED fit. However there are two complications with 
this approach: First, it depends on the assumptions that 
are made to obtain the mock images (dust modeling, stellar 
populations synthesis modeling). Second, the reduction of 
the mock images would have to be done in the same way 
as for the observed data, which is non-trivial. Therefore we 
have decided to conduct a simpler comparison within the 
limits of the quoted observational errors. 

The characteristic shape of the SHM relation is a 
consequence of several different physical processes, that 
prevent the gas in the dark matter halo from cooling 
(or reheat cold gas) and forming stars. The contribution 
of each process varies for systems of different mass. In 
low-mass systems, feedback from supernova-driven winds 
can be very effective as winds can easily expel large 
amou nts of gas from haloes with low escape velocities 
fe.g.|Larsodll974l;lDekel fe Silk|[l98ej ; lOppenheimer fc Pavel 



2008 ; IPuchwein fe SpringellboTa T Furthermore the amount 



of cooled gas is redu ced by the presence of a UV and 
an X-rav background (Navarro fc Steinmet j 1 19971 ; iGnedinl 
l200d ; iBenson et aHl2002l;lHambrick et alj|201lf) and by cos - 
mic rays l|Jubelgas et al.ll2008l ; IWadepuhl fc Springelll201l[ ). 
These effects have different efficiencies for different halo 
masses, such that the low-mass slope of the SHM relation 
can change at the halo mass, where the efficiencies of the 
two processes cross. However, with the current data it is 
very hard to find this transition mass, as the error estimates 
are even larger than the possible variations in the slope. 

In massive haloes the star formation efficiency may 
be lowered by feedback from active ga l actic nuclei (AGN 
Springel et ail 120051; ICroton et al.ll2006l; iBower et al-lboOrj; 



Somerville et al. |2008|" Hambrick et al.ll201ll; iTevssier etlnl 



201ll ; lMartizzi et al.ll2012l ; |Puchwein fe Sprmgel2012l ) which 
can heat the gas in the halo and prevent it from cool- 
ing. Another effect that may be important for massive 
galaxies is gravitationa l heating |Khochfar fc Ostrikerll200Si ; 
I Johansson et al.ll2009l ). In hydrodynamical zoom- in simu- 
lations this effect is naturally included and results in a 
trend of lower galaxy formation efficiencies for more mas- 
sive halo es, although still at levels which are sig nificantly 
too high (|Oser et alj|20ld : iHirschmann et alj|2012h . In this 
framework, the characteristic halo mass for which galaxy for- 
mation is most efficient is the mass where the combination 
of the various physical processes that reduce the efficiency 



is minimal. For massive galaxies the amount of in-situ stars 
is affected largely by processes that are effective for more 
massive systems (AGN feedback and gravitational heating), 
while the amount of ex-situ stars is strongly affected by pro- 
cesses that are effective for systems of lower mass (such as 
SN feedback and an ionizing background radiation). 

Although the only observational constraints of our 
model are the observed SMFs from z — to 4, our predic- 
tions for the evolution of the cosmic SFR density are in very 
good agreement with direct observation. This indicates that 
the evolution of the SMFs is consistent with the evolution 
of the cosmic SFR density, although the integrated observed 
SFR density was previously foun d to be higher than the ob - 
served stellar mass density (e.g. I Hopkins fc Beacoml [200^ 1. 
There are two effects that help to reconcile this tension. 
First, we employ the recent measurements by S12, who find 
a considerably steeper low-mass slope for the SMFs. This 
results in a larger amount of stellar mass in the low-mass 
end and thus increases the observed stellar mass density. 
Second, we include the effects of observational mass error 
in our model. As a consequence we find an evolution in the 
high-mass end of the 'true' SMFs, especially at low red- 
shift. When both effects are taken into account, the tension 
between the evolution of the stellar mass density and the 
integrated SFR density is considerably lowered. 



6 SUMMARY 

In this paper, we present a galaxy evolution model that as- 
sumes a tight relation between the stellar mass of a galaxy 
and the mass of the dark matter halo in which it resides. 
With this approach we populate haloes with galaxies over 
the entire cosmic period from z ~ 4 to the present, using the 
TV-body Millennium simulations and observed SMFs. Unlike 
the classical abundance matching technique our new method 
is self-consistent, taking into account that subhaloes are ac- 
creted at earlier times and with higher masses than when 
they are observed. We employ a redshift-dependent param- 
eterization of the stellar-to-halo mass relation, and popu- 
late haloes and subhaloes with galaxies requiring that the 
observed SMFs at different redshifts be reproduced simulta- 
neously. 

We show that when observational mass errors are not 
taken into account, no consistent stellar mass evolution of 
massive galaxies can be obtained, as the central galaxy of 
a typical halo of mass M V1I = 1O 14 M0 would lose ~ 25 per 
cent of its stellar mass from z = 1 to 0. When we include 
observational mass errors, this tension vanishes and massive 
galaxies also have a physically consistent stellar mass evo- 
lution. We find that the halo mass where galaxy formation 
is most efficient is log(M v i r /M0) ~ 11.8 at z = and in- 
creases with increasing redshift to log(M v i r /Af0) ~ 12.5 at 
z — 4. Thus at high redshift, the peak efficiency of convert- 
ing baryons into stars occurs in more massive haloes. The 
normalization of the SHM ratio decreases with increasing 
redshift: While at z = the most efficient haloes have con- 
verted ~ 23 per cent of their baryons into stars, at z — 4 
the peak efficiency is less than ~ 10 per cent. 

We further compare our results to those of hydrody- 
namical zoom-in simulations and semi-analytic models and 
find that most of the newer models show good agreement 
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with the efficiencies needed to reproduce the local observed 
SMFs. However, almost all have efficiencies which are far 
too high at earlier epochs, indicating that they form their 
stars too early. At late times star formation is quenched, so 
that most of the stellar mass today is already present at high 
redshift. In the semi-analytic models the normalization and 
the slopes of the SHM relation are in good agreement with 
ours, but the halo mass where the efficiency peaks moves to 
lower rather than to higher masses with increasing redshift. 

We then combine our relation with merger trees ex- 
tracted from the simulations in order to predict the mean 
assembly histories of the stellar mass components within 
dark matter haloes (in situ and ex situ formation). We av- 
erage over a sample of haloes with the same mass at z — 0, 
such that our model yields the average stellar mass growth 
of central galaxies as a function of halo mass. We find that 
in low-mass systems the halo grows much more slowly than 
its central galaxy. For a system with M V - 1T = 10 12 M Q , half 
of the stellar mass of the central galaxy forms after redshift 
z = 0.7 while half of the virial mass was already assembled 
by 2 = 1.3. For massive systems the trends are opposite at 
low redshift: the fractional growth of haloes is much stronger 
than stellar growth, while a large fraction of the stellar mass 
of massive galaxies is already in place early. For a system 
with M v i r = 1O 12 M , half of the halo has assembled only 
by z = 1.0, but half of the stellar mass was already in place 
by 2 = 1.8. In very massive haloes with M V - 1T = 10 12 Mq the 
galaxies start to form very early with 25 per cent of their 
stellar mass in place already by z — 3.0. However, it takes 
them a long time to fully assemble: Half of their stellar mass 
is in place only by z = 1.6, while 75 per cent of the stellar 
mass is in place by z = 0.3. 

We then go further and identify the progenitors of a 
given halo at earlier epochs and derive the amount of stel- 
lar mass that has been accreted onto the central galaxy. 
Employing a recipe to account for mass loss from individual 
stars, we then compute the stellar mass that has been formed 
through star formation and convert this into a SFR. In this 
way we derive the SFR history as a function of halo mass. 
We find that the star formation history of central galaxies 
peaks at a certain redshift which depends on halo mass: for 
massive systems the peak occurs at high redshift while for 
low-mass systems it has not yet been reached by z — 0. Af- 
ter the peak has been reached, the SFR quickly drops, such 
that massive galaxies are quenched at z ~ 2. On the other 
hand, the stellar accretion rates increase with time for all 
systems. Massive systems have a high stellar accretion rate 
that dominates the total growth rate at late times, while for 
low-mass systems, accretion is insignificant at all times. The 
halo mass at which the accretion rate and the SFR are equal 
at z = is log(M vir /M Q ) ~ 12.7. 

Finally we demonstrate that our inferred SFR histories 
are in agreement with observation. We compute the SFR for 
all systems in the simulation and derive the evolution of the 
cosmic SFR density. Moreover we compute the specific SFR 
history for two different stellar mass ranges. We find that our 
choice for the parameterization of the redshift-dependent 
SHM relation is consistent with observational constraints, 
and that our inferred recipe for the mass dependent star 
formation histories gives reliable estimates. 
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